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Abstract — Applying the Riemann geometric maciiinery of 
volume estimates in terms of curvature, bounds for the minimal 
distance of packings/codes in the Grassmann and Stiefel mani- 
folds will be derived and analyzed. In the context of space time 
block codes this leads to a monotonically increasing minimal 
distance lower bound as a function of the block length. This 
advocates large block lengths for the code design. 

Index Terms — Sphere packings, space-time codes, Gilbert- 
Varshamov/Hamming bounds, Stiefel/Grassmann manifold 



I. Introduction 

This work is inspired by Barg and Nogin's paper [1] for 
asymptotic packing bounds in the Grassmann manifold, based 
on an asymptotic expression for the volume of metric balls. 
The basic estimates defining the bounds are given by the well 
known Gilbert- Varshamov and Hamming (or sphere packing) 
inequalities: In a compact manifold M without boundary 
furnished with a topological metric d, let us denote the volume 
of the metric ball of radius S as vol Bd{5) (this quantity is 
presupposed to be independent of its center). Then for any 
given do there exists a packing (or code) C C M with the 
prescribed minimal distance do and cardinality \C\ such that 



vol M 



< vo\Bd{do) (Gilbert- Varshamov) 



while for any packing/code C C M with data (do, \C\) 



(1) 



volBd(ido) < 



vol M 



(Hamming) 



(2) 



holds. 

Taking for M the complex Grassmann manifold Gj^^ of 
k dimensional complex subspaces of C", Barg and Nogin 
derived closed form expressions 



vol BdiS) 




2nA:+o(n) 



2n/c+o(n) 



(geodesic distance) 
(chordal distance ) 



(3) 
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as n — > oo, leading to 
/ 1 



arcsm 



/ 1 



2fi/fc 



< do < 



2fc 1 - 1 - 



(4) 



(5) 



for geodesic, respectively chordal distance (defined later on), 
whereas R denotes the rate 



i?=^log2|C| 



(6) 



Furthermore Han and Rosenthal [2] recently derived upper 
bounds on the minimal distance (more general: on the diversity 
of space time codes) for packings on the unitary group U(n). 

A general capacity and performance analysis of space time 
codes in Rayleigh flat fading MIMO scenarios without channel 
state information at the transmitter [3], [4], [5], [6] revealed 
that the appropriate coding spaces are indeed 

• the (scaled) complex Grassmann manifold Gj?„ (set of 
fc dimensional linear subspaces of C"), if the channel is 
unknown at the receiver 

• the (scaled) complex Stiefel manifold (set of k 
orthonormal vectors in C") if the channel is known at 
the receiver 

Here k corresponds to the number of transmit antennas and n 
to the block length of the codes and the work in [1] refers to 



G 



k,n 



as 71 



oo while [2] refers to Vc^ as k 



The aim of this work is to close the gap between those 
two results by deriving bounds on the minimal distance for 
codes/packings in Gj?„, V;f„ for arbritrary (fc, n) (section III): 
Applying the bounds (1), (2) with equality, the main task is 
to solve the equation 



Bd{S) 



(7) 



for (minimal) distances S in G^„, V^'^^, with respect to 
some appropriate distance measure d. To this end volume 
estimates for the volume of (small) balls Bd{5) induced by 
curvature bounds for Gj? „ and V^^ come into play. Associated 
comparison spaces with constant curvature and simple volume 
forms provide bounds for Bd{S). In particular the lower bound 
turns out to permit a simple closed form expression with 
respect to (fc,n). Its analysis culminates in Theorem IV. 1 for 
the geodesic minimal distance lower bound and Corollary IV.2 
for the minimal distance do of the corresponding space time 
codes. Surprisingly it turns out, that the minimal distance do 



2 



grows at least proportional to ^/n, while keeping the rate and and 
the transmit power per time step constant. That is, increasing 
the block length enhances the possible minimal distance, thus 
in coding spaces with large block lengths there exists codes 
with potentially better error performance than in 'small' coding 
spaces. Since most of the space time coding research efforts in 
the literature deal with small dimensional coding spaces such 
as C/(fc) (e.g. [6]), future research in the more general „, 
^kn promises performance gains. 

Apart from space time codes recent developments in the 
design of space frequency codes [7], [8] also indicate that 
the relevant coding spaces are subspaces of large dimensional 
Stiefel and Grassmann manifolds. Thus the achieved results 
here may be of considerable importance for space frequency 
code design. and 

This article proceeds as follows. Section 11 deals with nota- 
tional conventions and basic definitions concerning the Stiefel 
and Grassmann manifolds (the coding spaces for space time or 
space frequency codes). In section III explicit bounds for the 
minimal distance will be calculated and compared to results 
obtained elsewhere. Further analysis on the lower bound will 
be performed in section IV, culminating in Theorem IV. 1 . Its 
implications for the minimal distance in coding theory will 
be pointed out in Corollary IV.2. Finally section V gives a 
summary of the results. 

II. The complex Stiefel and Grassmann manifolds 

The complex Stiefel and Grassmann manifolds together 
with their topological metrics (coding distance function in the 
language of coding theory) considered in this work constitute 
the focus of this section. For the analysis in later sections we 
also need some explicit curvature computations and rigorous 
proofs, which can be found in the appendices A and B. 

Readers who are mainly interested in the results concerning 
packings/coding and who are willing to accept the (quite 
standard) differential geometric facts can read this section 
without reference to the appendices, where further details can 
be found. 

A survey of the geometry of the real Stiefel and Grassmann 
manifolds aimed at non- specialists can be found in [9]^' holds 
and for an elementary introduction to differential geometric 
concepts see e.g. [10]. 

A. The Stiefel manifold V^n 
The (complex) Stiefel manifold 

V/;„ := {$ e I = 1} (8) 



(1 denotes the identity matrix) can be equipped with the struc- 
ture of an [/(n)-normal homogeneous space, which justifies 
the coset representation 

V^„ = C/H/(Jc/(„%)) , $ = (9) 

($ e U{n)), in particular 

dime Vfc^„ = dime C/(n) - dime C/(n - fc) = fc(n - |) (10) 

^^The complex case considered here is similar to the real case but in some 
places certain peculiarities of the complex structure come into play 



vol yfe^„ = vol U{n)/ vol U{n - k) 



= n 

i=n—k-\-l 



n 



27r* 



(11) 



i=n—k-\-l 



{^-ly. 



For V^j-^ as a Riemannian manifold the concept of geodesies 
and geodesic distance can be apphed to obtain a canonical 
distance measure r^: Denoting the tangent space of the unitary 
group U{n) by u(n) consisting of skew-Hermitian n-by-n 
matrices, tangents of V^^^ may be represented as 

u(n) bX=(^^ , A e u(fc), B e C^"-'^)^'^ (12) 



V\2 



^II^L' = ^II^L' + l|S||p' (13) 



is the squared geodesic length of the geodesic connecting 
* = (S) e F«„ with $ = (cxpX)(i) G T/«„. Here 
exp denotes the matrix exponential and the geodesic distance 
between arbitrary points e V^^ follows from the 

isometric transformation ^ = and $ = 

The canonical embedding (8) of V^^ into the vector space 
^£nxk^ <•, •>e) motivates the definition of another topologi- 
cal ('chordal') metric/distance 



:= 11$ - , * e V^,^ 



(14) 



which is important for space time coding, where it represents 
the decision criterion at the maximum-hkehhood receiver, if 
the channel is known at the receiver ('coherent' case), see 
[4], [6]. Note, that d^ is entirely different from the geodesic 
distance r^. Nevertheless we have^' 

Proposition II.l 

For k = n or k < ^ the metrics dX and are locally 
equivalent, thus in sufficiently small neighborhoods there 
exist constants > 0, > such that 



(15) 



This equivalence Unks the abstract (geodesic) sphere packing 
problem to space time coding. The restriction to the cases 
k = n and fc < J is mainly for convenience, since the main 
analysis will concentrate on A; ^ n. 

Proof: Lemma B.l, B. 2, B. 3 in Appendix B ■ 

Remark n.2 

While it is an easy exercise to find = (Lemma B.l), 

no concrete values for have been obtained rigorously. 
However, for /c = 2, n = 4, 6, 8 numerical simulations led 



to a 



2-0.9- 



^'At first sight the proposition seems obvious, but one has to take into 
account that dr is expressed in terms of * 6 ^kn' whUe is expressed 
in terms of the space of tangents and these two spaces are linked by the matrix 
exponential which can not be written in closed form compare Appendix B. 
Furthermore unlike d^, is NOT induced by (geodesies with respect to) 
the seemingly canonical embedding C C"^*', compare Appendix A-I.2 
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B. The Grassmann manifold „ 
The (complex) Grassmann manifold 

Gin ■■=m\^^vpj 



of all fc-dimensional linear subspaces ($) of C" also carries 
the structure of a Z7(n)-normal homogeneous space with coset 
representation 

GL^c/H/(Ta(„%) . {^) = ^^'-' (17) 

(with $1 := (1,0)$) and 

dime = k{n - k) (18) 
The total volume of „ is 



vol G«„=voiy«„/vol U{k) 



n 



i=n— fe+1 



" 27r^- (19) 

11 (n - 1M 



Tangents become 



with squared geodesic length = HiJUp, but there is an 

alternative notation in terms of the vector of principal angles 1? 
between subspaces: To simplify matters let us assume k < n/2 
whenever we are in contact with the Grassmann manifold. This 
is no restriction, since for /c > n/2 we can always switch to the 
orthogonal complement. Then there are precisely k principal 
angles between the subspaces ((J)) and ((cxpX)(J)). 
Performing a singular value decomposition on the tangents 
(20) one obtains (compare A-I.3). ||B||p = ||i?||i,2 thus the 
geodesic distance r*^ between ((J)) and ((expX) (J)) reads 



\B\\ 



(21) 



As for V^^ there is also a different distance measure d'-^ in 
^ induced by the maximum-likelihood receiver, which can 
be derived from the following geometric picture: 

Spherical embedding: Unlike for the Stiefel manifold, there 
is no canonical embedding of Cj^^ into Euclidean space 
unless choosing a representing unitary frame $0 in each 
subspace ($) G G^„. Nevertheless there exists an interesting 
embedding of G^?^ into Euclidean space given in [11]: For 



projection 



there is an well-defined associated orthogonal 

$$t . — , (22) 



* : = 



of norm \\P — k/nl\\l = k(n — k)ln and tiP = k, which 
justifies the embedding 



'■{^k{n-k)/n) C R" 

This motivates the 'chordal' topological metric 

, (*)) := ||sint?|U2 = -1=\\P^ - P^IIp 



(23) 



(24) 



(4>, 4" € V^j^). Comparing rf"^ with the geodesic distance r*^ 
(21) between two subspaces we observe 



(16) Proposition n.3 



(25) 



whereas fi'^ = 1 and a'^ = f . 



III. Bounds for the minimal distance 

Now let us specialize the general packing/coding bounds 
(1),(2). Set 



(M,rf) := 



and (compare (10), (18)) 



D := dimia M 



k{2n-k), M = y«„ 



2k{n-k) 



(26) 



(27) 



for the two cases of interest. In the sequel other symbols Uke 
a are used genericaUy to denote or a'^ when specialized 
to the corresponding spaces V^^, G^„. Denote by 



v{r) := volB(r) 



(28) 



the volume of the geodesic ball of radius r in M, which is 
independent of its center by left invariance of the Riemannian 
metric. With this notation the Gilbert- Varshamov (1) and 
Hamming bound (2) for packings C on M can be compactiy 
rewritten as 

1 /vol M\ 1 /vol M^ 



1 /vol M\ „ 1 /vol M\ ^ 

or relaxed w.r.t. coding distances do using (15), (25) 

1 -i/volMx , 2 i/volM\ - , 

^0 := -^~\-^) <do< J^v-\^;nr) =■■ do (29b) 



So packing bounds are related to the coding bounds by simply 
setting a = = 1, thus replacing the (topological) metric 
distances by geodesic distances. Due to the rather difficult 
to obtain exphcit value for in (15) we focus on the 
packing bounds (29a) for most of the remaining analysis, 
keeping in mind the simple relationship between statements 
about packings and statements about space time coding. 

To obtain the desired bounds for the minimal distance 
provided by (29) we need closed form expressions for the 
volume V of small balls in M. As has been already indicated 
in the introduction, this is a difficult task in general: The 
canonical volume forms on G^„ and V^\^ are elaborate to 
calculate. Alternatively a common tool to compute volumes 
in Riemannian geometry arises from curvature, using Jacobi 
vector fields (see e.g. [12] for details). Unfortunately a direct 
apphcation can not be performed since we would have needed 
a diagonaUzation of XY — YX for each horizontal (compare 
A-1.1) = = 1 in u(r7.) written in closed form. But 
there are simple volume estimates which will be presented in 
III-A. In III-B the results will be compared to those already 
obtained in [1], [2], in a few (computational simple) cases. 
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A. Bishop/Gunther volume bounds 

The method for volume computations in Riemannian man- 
ifolds using Jacobi vector fields can be looked up in [12, 
theorem 3.101]. For k e E, let 

D-l 



(30) 



denote the volume of the geodesic ball of radius r in the 
manifold of constant curvature k and let k < k < k be defined 
by (compare (A.5), (A.4)) 

1 



K := 



min Ric(ei, e^) 



k:= max K(X,Y) 
11:^11=11^11=1 



(31) 



then we obtain monotone volume bounds Vi{r) < v{r) < 
Vu{r) for arbitrary < r < by 

vi{r) -.^v^ir) , vuir) ■.= v^{r) (32) 

From X < K ^ ■u''(r) < v^{r) and K > in M (A.6) 
we can further relax k to zero, which yields the simple upper 
volume bound 



(33) 



A lower volume bound comes from an upper bound k for K. 
Inserting tangents X, Y (12), (resp. (20)) into (A.4) subject to 
||X|| = ||F|| = 1 yields 



K{X, Y)<K = < 



2 (n?„,fc<n) 



.4 (G^,J 

Plugging this bounds into (29a) we we end up with 



2 iUik) = V^.) 



(34) 



.o = (^V(^)<ro<2K)-(^)^.o(35) 
With these settings an expUcit (Maple-) calculation revealed 



k\n 


k 


2k 3k 


Ak 




1 


[1.57,3.14] 


[1.06,3.49] [0.941,-1] [0.886,-1] 
{0.500, 1.05}{0.595, 1.40}{0.630, 1.71} 


2 


[1.58, -1] 


[1.38,-1] [1.32,- 
{0.771,-1} {0.909, 


-1] [1.29,- 
-1} {0.973, 


-1] 
-1} 


3 


[1.74,-1] 


[1.66,-1] [1.63,- 
{0.977,-1} {1.15,- 


-1] [1.61,- 
-1} {1.24,- 


-1| 
-1} 


4 


[1.92, -IJ 


[1.92, -IJ [1.89,- 
{1.15,-1} {1.36,- 


-IJ [1.88,- 
-1} {1.46,- 


-IJ 
-1} 



[ro, ro] for V^^^ and {ro, ro} for G^^ with respect to (35) for ii = 1 



fc\n 


k 


2k 3k 4k 


1 


[0.0031,0.0061] [0.0165,0.0330] [0.0223,0.0446] [0.0251,0.0502] 
{0.001, 0.002} {0.0055, 0.0110}{0.0098, 0.0197} 


2 


[0.0700,0.140] 


[0.172,0.348] [0.203,0.412] [0.217,0.441] 
{0.0341, 0.0682} {0.0877, 0.176} {0.122, 0.245} 


3 


[0.217,0.440] 


[0.416, 0.898J [0.467,1.04] [0.490,1.11] 
{0.122,0.246} {0.242,0.504} {0.309,0.664} 


4 


[0.403, 0.847J 


[0.678, -IJ [0.743, -IJ [0.771, -IJ 
{0.242,0.504} {0.422,0.992} {0.517,1.75} 



[r;o,roJ for V^^ and {ro,ro} for ^ with respect to (35) for il = 10 

(— 1 in the tables means, that Maple could not find a 
solution, due to approximation error/too large sphere radii). 



Observe that in the lower rate regime ro still grows with n 
in C^j^, but slowly decreases in V^^, while in the high rate 
regime ro is strict monotone with respect to n, expecting the 
intervals to become disjoint. 

So, while the high rate requirement is too restrictive, the 
results for low rates are unsatisfactory in part. But the general 
analysis of the lower bound in section IV will come up with 
interesting results, supporting this approach. To clarify the 
presentation let us sunmiarize the results so far in the 

Proposition III.l 

The inequalities (35) provide approximate bounds on 
the (geodesic) minimal distance for packings/space time 
codes on the Stiefel (coherent case) and Grassmann (non- 
coherent case) manifolds for any admissible {k, n). In par- 
ticular the lower bound in (35) is computational simple and 
guarantees the existence of corresponding packings/codes. 

Especially for the Stiefel manifold these explicitly calculated 
bounds appear to be new in the context of coherent space time 
coding. 

B. Comparisons with related results in the literature 

Han and Rosenthal [2] obtained bounds on the scaled 
chordal distance A := in the unitary case k = n, 
^k^k =U{k). Based on a numerically calculated exact volume 
they extracted three upper bounds. The following table shows 
their (best) upper bounds (2nd row) for A = d"^ in 
U{2) (in part relying on the results in [13]) for different rates 
i? (1st row) together with the upper bounds obtained here (3rd 
row) 



R 



2.29 2.79 



3.0 



3.16 



3.32 



3.45 



3.50 



4.98 



A I 0.675 0.619 0.597 0.580 0.558 0.542 0.535 0.327 
^ I 1.40 1.01 0.909 0.843 0.785 0.742 0.727 0.409 

Note, that equality in the (rough) estimate 7^7= < ^ 
(Lemma B.l) has been forced in the third row of the table to 
convert the geodesic distances computed by a Maple program 
into chordal distances. Consequently the bounds of [2] are 
tighter than the bounds obtained here, since in the case of 
unitary matrices there are more specialized (but less general) 
methods available to obtain bounds. 

As already stated in the introduction another (asymptotic) 
result has been obtained by Barg and Nogin. For the non- 
asymptotic case they presented an exact volume formula [1, 
eq. (11)] for regions in the (real and complex) Grassmann 
manifold. 



volB(r)=2'=|G«jn 



{n — i)\ 



0<i?i<---<l?fc<7r/2 
Il<'ll2<r 



ii[(i_l)!]2(„_fc_i)r 
d^i . ..d'dk'x 



(36) 



]J(sint?i)2("-2'=)+i n(sin2 j9j - sin^ ^if 

i=l j<l 
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which can be computed in polar coordinates, compare A-II. 
Although (36) is exact, it does not provide a closed form for 
varying dimensions. Moreover the computations are elaborate 
compared with the ones done here, such that the evaluation of 
(29) become intractable. 

IV. Analysis of the lower bound 

The lower bound for the (geodesic) minimal distance tq 
guarantees the existence of corresponding packings/codes. Due 
to (33) we can expUcitly solve the lower bound in (35): 

D„,k _ J_ vol M 



(37) 



with -D„,fc defined as D in (27), thus 

f 1 

12, G«„ 

Then 



(27') 



Theorem IV.l 

The (geodesic) minimal distance tq in M can be lower 
bounded by 

to > (l) (38) 



with the right hand side monotonically increasing as a 
function of n for n> k. Asymptotically 



lim Tg 



2R/k 



(39) 



holds. 



In particular this establishes a monotonically increasing lower 
estimate for ro common for V^^ and G^^, which is not 
obvious from the picture drawn from the explicit calculations 
of ro for rate R = 1 in the previous sections. Of course, 
the theorem also holds for the (topological) minimal distance 
do = ^ro, connecting this result with space time coding 
theory. 



^)'/''-Mhen ro = a-6 



Proof: 

Set a := 2-''^/^'--" and b := (- ^ 
and from (27') a is monotonically increasing as a function of 
n with lim„^(x) a = 2"^/^'=. 
For b we show 

b > 1 (40) 



and 



lim b = Vk 



(41) 



and the theorem follows. 

For the two cases of interest 6^" '= is given as (using 
(A.10),(ll),(19)) 

vol F«„ ^ r(fc(2n - k)/2 + 1) 



|^fc(2r!.-fc) I 

vol 

|_g2fc(n-fc) I 



n 



i—n — k-\- 



r(2fc(n-fc)/2 + i)nLr(j) 



n 



i—n—k-\- 



for n > k, resp. n > k + 1. The proof of (40) relies on the 
simple estimate 

^^^i±^ = {m + l){m + 2)- ... -M 
r(m+l) ^ ^ (43) 

for m,M G ilN,M - m e M. Since Dn,k > it suffices to 
show Bk,n := b^"-'' > 1 for all admissible (fc,n). This will 
be proven by induction over k and n. 

1) = 20Fr(3/2) = TT > 1 

2) Induction over k 

M 



Bk+i,k+i _ 2^r(^'/2 + + 1/2) +1) 



Bk,k 



r(fc2/2 + i)r(fc + i) 

(fcV2 + l)'=+V2 
k\ 



> 
(43) 

> 2y^y/k^/2 + 1 

fe2/2+l>fe 

> VQtt> 1 

k=l 

Thus Bk,k > 1 Vfe>i 
3) Induction over n> k 

m 

Bk,n+i T{k/2{2n-k) + k + l)T{n-k+l) 



Bk,n 



T{k/2{2n-k)+l) T{n+1) 



_ (m+l)- ... ■ (m + k) 
(43) (m + 1) • ... ■ (m + k) 



> 1 



since m - m = n{k - 1) - P/2 + k > k'^/2 > 0. 

n>k 

Thus for every fc > 1 we have Bk,n > 1 V„>fe. 



_ r(2)r(i) _ 
'i>2 - r(2) - 
2) Induction over k 



1) Bi 



Bk+i,k+2 _ nk+2) .^ n'^2 _ 
Bk^k+i T{k+i) ^ ^ ^n-i'rw 

Thus Bk,k+i = 1 Vfc>i 
3) Induction over n>k + l 

rh 

Bk,n+i ^ T{k{n-k) + k + l)T{n-k+l) 
Bk,n ~ r(fc(n-_fc)+l) r(n + l) 

m 

(m + 1) • . . . • (m + fc) ^ ^ 
(43) (m + 1) • ... • (m + fc) ~ 

since m — m = n{k — 1) - k"^ + k > k — 1 > 

n>k+l 

and it foUows for every fc > 1, that Bk^n ^ 1 ^n>fc+i 
as desired. 

Let us now prove (41). At first (20r)''/^"''= — > 1 and 

n— *-oo 

nLi rOy/^" '' — > 1 holds. So it remains the evaluation 



(42) 
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of 



Stirling's formula reads either ('~' denotes asymptotic equiv- 
alence) 



T{m+l)r~.^/2^i^—j or r(m) 



/27r /my 
m V e / 



and by Dn,k ^ 2nfc we deduce 



r(D„,fc/2 + 1) X 1/^..-= ( \/27r£'„,fc/2y^ 



(nr=„-.+ir(i) 



1 li=n- 



fe+1 



27r / i 



(27r)* 



\/2e 
1 



1 



In ■ 



— W2ek 

This proves (41) ■ 

A. Final remarks and application to coding theory 

A remarkable coincidence arises from Barg/Nogin's results 
for the chordal distance in G^^. Denoting the lower bound in 
(5) by ^0 we find 

^0 = lim ro (44) 

n — *oo 

therefore, the geodesic lower bound Tq obtained from the flat 
geodesic volume estimate v^{r) < v{r) asymptotically equals 
the (asymptotic) exact chordal lower bound (5). This seems 
reasonable since in flat space the geodesic distance coincides 
with the EucUdean (chordal) distance. 

Apart from the asymptotics, let us consider the lower bound 
(38) of Theorem IV. 1. It guarantees the existence of pack- 
ings with minimal distance ro bounded monotonically from 
below in V^\^, resp. G^^, when n grows. In coding theory 



^k,n — 



k ^k,n 1 " 



k 'k^n . " resp 
represent the coding spaces for space time block codes for 
the Rayleigh flat fading channel unknown to the transmitter 
and known, resp. unknown, channel at the receiver. The factor 
serves as a constraint, holding the transmit power at 
each time step constant for different choices of (fc,n), thus 
provide a fair comparison of codes from different coding 
spaces. In a Riemannian manifold M with metric g the 
mapping (AM, g) i — > (M, X'^g), A > 0, is isometric, leading 
innmediately to the scaled geodesic minimal distance fo = Xtq. 



With respect to the coding distances d we obtain instead 



(45) 



whereas (p > 1 denoting the signal to noise ratio) 

is (a lower bound of) the first order term (the so called 
diversity sum, our metric here) in the expansion of the Chernov 
bound for the pairwise error probability, compare [4, formulas 
(17)(18)(19)(20)]: The factor of i for stems from the 

slightly different 'effective' transmit power g := ^i^^J^/i^^ 
compared to the known channel effective "transmit power 
g := H, satisfying ^g < g < g, whereas p > l,n > k 
is understood. Collecting aU formulas we finaUy infer from 
Theorem IV. 1 : 

Corollary TV.2 

Given p > 1 and n > 2k, there exist space time block 
codes with minimal distance do lower bounded by 



(46) 



whereas o. is d€tcvTfiiTicd by (15), resp. (25), Dn 

is defined 

in (27) (resp. (27')), and ji in (45). 

Thus the performance (which scales with d^) potentially 

increases monotonically at least proportionally to ^. 

The last statement in the corollary follows from the observa- 
tion, that the diversity (essentially the inverse of the Chernov 
bound for the pairwise error probability) as a basic perfor- 
mance measure for space time codes [4] is a homogenous 
polynomial. The first order term coincides with the metric d^, 
while all higher order terms scale with a power of when 
code design is interpreted as a constrained packing problem 
(considering the higher order terms as constraints according 
to a normalized distance distribution). 



V. Conclusions 

The framework in [1], [2] had been successfully generalized 
to flie Stiefel manifold V^^^, n > fc, and to oo » n > 

fc/2 using the completely different method of Riemannian 
volume bounds (Proposition 111.1). Unlike the exact volume 
formula the lower bound can be relatively simple analyzed as a 
function of (fc, n) for both and V^.^, leading to Theorem 
IV. 1, resp. Corollary IV.2. Although the used estimates were 
quite conservative they apply (in principle) in any Riemannian 
homogeneous spaces. 

The connection to the coding theory of space time block 
codes advocates further efforts in finding codes in the spaces 
V;f„, resp. Gj?„ for n much larger than k. Since the minimal 
distances grow proportionally to while the transmit power 
per time step remains constant, there is a considerable perfor- 
mance impact to expect, when coding in V^^ (resp. G^ „) as 



7 



opposed to coding in U{ky\ Furthermore, as already pointed 
out in the introduction the developments in space frequency 
coding indicate, that the relevant coding spaces are subsets in 
some Vjpyj (resp. G^^) whereas the number of subcarriers n 
satisfies n ^ k, thus the results proven here may apply to 
space frequency codes as well. 
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Appendix A 
Differential geometric calculations 

/. U {n)-normal homogeneous spaces 

For the theoretical background, common notation and cur- 
vature formulas I refer to [15] as a reference. 
1 ) The unitary group: 

C/(n) = {!> e C"""" I = 1} (A.l) 

is a compact, connected Lie group and a real manifold of 
dimension dim^ — r? . The corresponding Lie algebra (i.e. 
the tangent space of U{n) at 1) is 

u(n) = {X e C"""" I = -X} (A.2) 

and the matrix exponential exp maps u(n) into U (n). On u(n) 
the (bi-invariant) Riemannian metric for V{n) is defined as 

<X,y>=^trXtF (A.3) 

thus <X,X>= \ \X\^ (Frobenius norm) holds. 

A manifold M is a C/(n) -homogeneous space, if there is 
a transitive ?7(n) action on M such that M = U{n)/H for 
some isotropy subgroup H c U{n). If \) c u(n) denotes 
the Lie algebra of H there is a canonical decomposition of 
tangent vectors f} ® ()^ = u(n) 9 X = X" -|- X and we 
can identify tangents of M with so called 'horizontal' tangent 
vectors X e f)"^. With this identification M is called normal 
homogeneous. 

Then the sectional curvature K and the Ricci curvature Rio 
of M axe given as 

K{X,Y)^\\\\X,Y]r + \\\[X,Y]Y (A.4) 

Ric(ei, Ci) = ^ K{ei, ej) (A.5) 
j 

whereas [X, Y] = XY—YX, X and Y are normalized tangent 
vectors and {ci} denotes a orthonormal base in ()^. Note that 
the sectional curvature K is always non-negative 

K>0 (A.6) 

^'Note that this does not contradict the (converse) conclusions in [14], which 
do not apply here: The error probability computations done there with respect 
to increasing block length n — » oo are constrained by a fixed total number of 
code symbols sent. This is a different scenario, not relevant for the analysis 
performed here. 



2) Supplements for the Stiefel manifold: The (complex) 
Stiefel manifold (8) is canonically a f/(n)-normal homoge- 
neous space: The canonical left multiplication of A;-frames 
in C" by unitary n x n matrices transforms each pair of 
fc-frames into each other. Thus the group action of U{n) 
on Vj^^ is transitive with isotropy group = ( J [/(^^.j.) ) 
and establishes the canonical diffeomorphism (9). Then 1) = 
(o u(n-fe) ) and tangents X £ i)^ have the form (12), and (13) 
follows for the geodesic distance r^. Note that this distance 
is not induced by the length of the geodesies obtained from 
the canonical embedding of VJ?„ into C"^*^, compare [9] in 
the real case and additionally [15, Example 6.61(b)] in the 
complex case. 

3) Supplements for the Grassmann manifold: The (com- 
plex) Grassmann manifold (16) carries the structure of a 
C/(n)-normal homogeneous space by forgetting not only the 
orthogonal complement of $ e C/ (n) (which has been done 
for V^^) but also the particular choice of the spanning k- 

frame. Thus H = ^^^^ u{n-k) ) ^hi^ '^^^^ to (17). Note 

that the coordinate representation ($) = ^ holds only 
locally in general, but it turns out, that this representation 
covers all but a set of measure zero, hence we abandon this 
distinction between local and global properties in this work 
and drop the distinction between G^^ and its coordinate 
domain. Calculating i)^ leads to tangents of the form (20). 

Given two elements ($) , e G^? „ then the k sta- 
tionary angles < iJi < ■ ■ ■ < 'Ok < 7i'/2 between 
($) and {^) are defined successively by the critical values 
Bxccos\< Vi,Wi >|, i = l,...,k (in increasing order), of 
(Ujw) I — > arccos|<w, w>| where the unit vectors w,w vary 
over {vi, . . . , Vi-i}-^ C ($), respectively {wi, . . . , Wi-i}-^ C 
{9). It is well known that the stationary angles can be 
computed by the formula (any representing fc-frame will do) 

cosi?i = (Ti($%), i = l,...,A; (A.7) 

whereas (Tj(M), i = 1, . . . , fc denotes the i-th singular value 
of the matrix M in decreasing order. 

Given a tangent X = (°"o')' ^ ^ ([|(«-fc)xfc^ 
singular value decomposition B = VT,W'^ = ViSW\ 
V = iVi,V2) e U{n- k), W e U{k), E = (g), 
S = diag(o-i,...,crfe), yields X = UAW with U = 
(T^.T^J> A = (^[]), and ^ =^(J^"o'')- From this 
one calculates (cxp X) (k) = ( 1^^™''?'^+ h and cosi? = 

^ ' \Vi{sin S)W^ ) 

cr( ( 1 ) (expX) ( J ) ) = cos 5, thus i? = cr and (21) follows. 

The space of orthogonal projections Ily := {P$ | $ € 
^kn) (compare (22)) can be identified with G^^. In particular 
we have Ily = 11^ with 

Efe := jPe C"""" = P,P^ = P,trP = fc, 

o . (A.8) 

\P -klnVg = k(n-k)ln\ 

as one can see by picking an appropriate representative $ e 
($) (e.g. $ = ( J ) due to invariance of II/j under the left and 
right unitary action). Since each P g 11^ is Hermitian with 

2 1 

constant trace, 11^ is canonically a real submanifold of E," ~ , 
the constant norm justifies the embedding (23) 
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//. Volume computations 

1) Total volume: The unitary group U{n) C GL(n, C) can 
be equipped with the induced Lebesgue measure from the 
ambient space R^" . The Stiefel manifold inherits its volume 
measure from its total space U{n): We get from the familiar 
volume formulas 



: vol 5 



m—l 



r(m/2) 
:=vol B™= |5"-H/m 



(A.9) 
(A.10) 



for the unit sphere 5™^^ and the unit ball in E.™ and 
the canonical homogeneous family 5*^™^^ = U(m)/U{m — 
1) the following recursive formula vol [/(I) = \S^\ — 2tt, 
vol [/(to) = |52™-1|vo1[/(to - 1), and therefore (11) and 
(19). 

2) Volume for regions in G^^ •" The volume formula for 
regions in the complex Grassmann manifold will be derived, 
based on [1, eq. (11)]''\ 

Starting with formula [16, (A. 18)] for the distribution of 
eigenvalues Ai, with Ai = cos^ -di of ( i,o ) $<1>^ ( o ) we obtain 
the volume density as the marginal density 

k 

w = C{k, n) ■ - Xi)^-^' Y[{Xi - Xjf • dAi . . . dXk 

i=i j<i 



C{k,n) ■ 2''k\ ■ ]J(sin^?i)'^"""'^+'cos 
JJ(sin2^j -sin^ dif. 



3<l 



Ml... Mk 

(A. 11) 



whereas the Jacobi determinant 2*^ J^jLj^ sini?.; cosiJi of the 
mapping A i — > -d has been introduced in order to express the 
volume density in terms of ??, and k\ establishes the ordering 
condition on the (open) simplex 9 = {0<i?i <•••< 
'dk < 7i'/2} of stationary angles. The constant C is just a 
normalization factor, which reads in our case 



C{k,n) 



\GZJ 
k\ 



{n — i)l 



(A. 12) 



(without the factor [Gl'^l this would give the Haar measure 
used in [1] on „). The volume of sufficiently small geodesic 
balls is now given as 



volB(r) = / w 

ien{||i?||2<r} 

= fdpf 

Jo Ja 



w(i?(p,a,/3))|det J^(p,a,/3)| 



da dPi . . . dpk-2 



(A. 13) 



'•'Unfortunately (in their first paper version) their formula is not correct in 
the complex case (private conmiunication). Fortunately this does not affect the 
(asymptotic) results obtained in [1]. An erratum has already been produced, 
thus the derivation here is only for completeness of the presentation and the 
convenience of the reader 



whereas (p, a, (3) denote (k dimensional) polar coordinates 

'&i= p cos (3k-2 ■ ■ ■ COS (3\ cos a 
i?2 = P cos (3k-2 ■ ■ ■ cos /?! sin a 



(A. 14) 



dk-i = cosPk-2SmPk- 
•dk = sm(3k-2 



The factor ^ removes the ordering condition on the simplex, 
such that the domain of angle integration is the whole region 
[0, 7r/2]'^~^. Eventually, denotes the Jacobi matrix of the 
coordinate transformation (p, a, /3) i — > 

Appendix B 
The local equivalence of d and r in V^^ 

In this appendix the proof of Proposition II. 1 will be carried 
out. Let us recall, what we want to show. Given in 
the complex Stiefel manifold V^,^ C C"^*^, the topological 
distance d motivated from coding theory is given as d = 
11$ — ^'IIf (we drop the upper index 'V in this appendix). 

At the same time, locally there is a unique geodesic 7 in 
V^,-^ joining $ and ^P, and the geodesic distance r is simply 

defined as its length L = /o^||7(t)||FC?i, 7(<) being the parallel 
transported horizontal tangent vector X{'y{t)) along 7. Thus 
we obtain r = ||X(7(0))||f. Since both d and r are invariant 
under the action of the U{n) we can set 4* = (p) without 
loss of generality. Recalling the general form X = ^ , 

A e u(fc), B G £{n-k)xk^ horizontal tangent vectors in 
u(n) (12) we arrive at 

d' = \\'i>-ih)\\^ 



1, 



= ^wxw 



1, 



(B.l) 
(B.2) 



whereas $ = cxp X ( J ). Unlike the case ^ = (representing 
tangents for G^? ^) there is no closed form expression for $ in 
terms of X in general (compare [9]), so it remains a non-trivial 
task to find constants a,p > satisfying 



I3d<r<ad 



(15') 



expressing the equivalence of d and r. 

To begin with the easy cases, the constant /3 is easily found, 
as well as a when k = n: Both are simple consequences of 
the two sided inequality sin a; < x < (7r/2)sina;, whereas 
X e [0, 7r/2] is understood in the second inequality. 

Lemma B.l 

In 14*^^ ■^d< r always holds, thus we have P = 
Proof: 

Since X e u(n) there exist V G U{n) such that X = 
Vdiag{i^)Vl thus r2 = ^ = (a,...,^„) e R". 

Now we can estimate as follows 

= 11(1 -expX)(S)||2< 111 -expXIL^ 



= 111- 



exp(diag(zO)||? = l]|l 



= 2^(l-cosC,) = 4 
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(since sin^ x/2<x'^ /^) ■ 

Lemma B.2 

If k = n then r < -^^d holds, thus a = 

Proof: 

k = n implies B = 0, X = A and we can estimate 

d' = \\l- exp = 4^sin2 ^ > ^Jaf = 

3 

(since a;^/4 < (7r^/4) sin^ a;/2 for x € [— tt, tt]), whereas la = 
i{a\, . . . , o„) denotes the vector of eigenvalues of A e u(n). 

■ 

The non-trivial task is to obtain some a > 0, when k < n. 
The rest of this section deals with this job. Let us assume 
k < J since this is the relevant case for the analysis in this 
work (the case A; > f should be similar). 

Let X = y + Z with y = 0), and Z = ( -^^ ), then 
we can write 



$ = exp X = (exp Z){^°) 



(B.3) 



since this is merely a factorization of $ = $ ( □ ) into a certain 
projection onto Cj^^ and the remaining 'phase' in U{k) 9 v. 
The first factor exp Z can be calculated in closed form: B has 
a singular value decomposition B = diag(i9^)ii for some 
V &U{n- k), u e U{k) and i?^ := {dk, ■ ■ ■ denotes 
the vector of principal angles (in decreasing order) between 

diag(cos ^^"^ ) — diag(sin^^) 



((o)) (^)- Setting U^l = I diag(sin,?i) diag(cosi>i) 

V 1 

we arrive at exp Z = ( q ^ ) f/,<,j ( q y )^ So we have achieved 
a quite explicit representation of $. In particular the principal 
k X fc-submatrix (/> = ( i o ) $ ( ^ ) reads 

<j) = u diag(cos 'd^)v}v (B.4) 

Now we can start estimating: 

<^ = \\^-{l)\\l = , 2{k-mttTci>) 



K 



(B.5) 



> 2[k-lY.[mc{u^vu)jjf - i^cos^^,) 



Writing U{k) 3 v = exp A, A G u{k) with eigenvalues 
la = i{ai, . . . Qk) of A we have 9lt{u^vu)jj > whenever 
d G [— 7r/2,7r/2]'^. Demanding this mild locality restriction we 
get [me{uhu)jj]^ < me{u^vu)jj, thus 'Ej[^^iu^vu)jj]^ < 
y{ttv{u^vu) = fHetrw = ^ . cosaj and therefore 



(f>2Y^ sin 



2 "J 



-^sin^ 1} 



(B.6) 



3 J 

All what remains to do in order to compare d with r is to 
find the link between A and A, respectively Y = ( q q ) and 
y = (^g).By(B.3) 



expy 



V 
1 



exp(— Z) expX 



(B.7) 



holds, thus our 'missing Unk' is given by the Baker-Campbell- 
Hausdorff formula expressing W G u(n) given by exp W = 
expf/ expy, ([/, V) e u(n)p x u{n)^ by 



W=V + 



')Udt 



r + 1 



E 



V Pi! qi\ P.< q.\ ) ' 



(B.8) 



fJi,...,p,^>0 

giv.9r>0 

Vi=i..^Pi+gi>0 



Pl^ h Pr + 1 



whereas /(z) = ^ and adj/ : F i — > [U, V]=UV - VU 
(see [17] for that particular representation of the BCH formula 
(Dynkin's formula in their terminology)). The second part of 
(B.8) is nothing but the term- wise integrated Taylor series 
expansion of the integrand. Following [17] the domain of 
definition u(n)g is the region of u(n) in which the tangent 
map of exp is regular. It is the complement of {?7 € 
u(n) I dct(adc/ -27r«Z'l) 0, = Z \ {0}} in u(n). In 
particular, u(n)g contains a connected neighborhood 



D{5q) = {Ue u(n) I \\U\\, < So} 



(B.9) 



of 0. Speciahzing to W = Y, U = -Z, V = X = Y + Z 
yields in multi-index notation (thus \p\ = J2iPi' = Hi ft') 



E 



(ad|^ ° ad«^ o . . . ° adP" ° ad«,-) (Z) 



P=(pi,...,Pr)>0 

g=(9i,...,'/r)>o 

p+<i>0 



(bl + 



(B.IO) 



Note that every term contributes at least some factor involving 
A (since &Az{Z) = 0), hence in the norm estimate 



||ad|^= ad«^o ...o adg-o adJ(Z)||p 
<2lpl+kl||2'||lfl+i||X||l«l 



(B.ll) 



X=Y+Z 



< 2H+M||Z||I''I+^^ I'^l ||y||j||Z 



ikl- 



the term corresponding to ? = has no counterpart in (B.IO) 
(resp. it is zero in (B.IO) already), therefore with \\X\\^ < 
5 (thus ||y||F, II-Z'IIf < 5) we can factor out one ||y||F and 
estimate 




)\\yh (B.12) 

<5i9i-i||y||p< 2i«i^i«i-i||y||F 
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and the k x k principal submatrix of (B.IO) of our interest 
satisfies 



A = A + C (B.13) 
\\C%<k\\A\\, (B.14) 

^ {\p\ + l)p\q\ ^^-^^^ 

P={pi,...,Pr)>0 

9=(Ql,...,q,)>0 
p+q>0 

It is possible to rewrite (B.15) such that we can prove the 

convergence of the multi-series, that is existence of k. given 
some sub-multi-indices pj, qj corresponding to some J C 



OD -J 



{l,...,r}letus setAj := 



and 



{ 4S)'" -' 



then 



(B.15) equals "fi'^r- with (J' denotes the set {1, . . . , r} \ 
J) 

r 

s=0 Jc{l,...,r} 

\J\=s 

I pj>l g>l p>l qj>l ) 



pj,=0 



qj'^ 

qj,=0 



(B.16) 

Now we can perform a rather rough estimate on the sums. We 
haveE,,>iM,. = (e^'-ir^ and Ep,,>i A,,, < (e^^-l)!''!, 
therefore (note that the sums in the brackets in (B.16) do not 
depend on the particular choice of J C {1, . . . , r} but only 
on its cardinality | J| = s) 



= [(e4* _ l)e2Y + [(e^^ - 1)6^^ 

(B.17) 

It is obvious, that we can choose a S < 6q sufficiently small, 
such that Kr < (^j.^i^t for any given t > 0, ensuring the 
convergence of (B.15). Setting in particular t = 1 yields k < 
^ — 1 < 1 and we obtain from (B.13) 

||A||p> (1-k)||A||p (B.18) 

Now we can proceed further with (B.6): 



and we have proven our final lenmia: 



Lemma B.3 

If k < ^ there exists a S < Sq, such that (1 — k) > 0, 
whereas Sq and 6 are determined by (B.9), resp. (B.17) 
demanding Kr < (r+i) - Then locally for r = \\X\\p < 5 
the relation r < 2{i-k.) ^ holds, thus a = 2(v-k) - 

This lemma fills the gap in formula (15). Of course, (1 — k) ~ 
1 would be optimal in this situation (observe the loss compared 



to a in Lemma B.2), which can be achieved by setting 
S <^ 1, with K decreasing the smaller S has been chosen. 
Unfortunately, the smaller we choose 6, the higher the required 
corresponding rate R ensuring the validity of Lemma B.3. For 
example, to obtain a numerical value of i?o ~ 1-4 (by formula 
(35) as a lower bound for the corresponding rate, with 6 — ro), 
which is still achievable for coding purposes in a practical 
setting, one needs values of (5 « 1.25, which is quite large 
in order to apply Lemma B.3, thus the estimates done here 
are far to rough to accomplish that. The importance of the 
lemma actually hes in the fact, that it proves the existence of 
some a > in (15). However, numerical simulations indicate 
that the real world behaves much better than the estimates. 
The histograms in Fig. 1 display 1 — k drawn from 1000 
random samples in ¥2^^, n = 4, 6, 8 for 5 = 1.25: Although 
there seems to be no rigorous and essentially sharper estimate 
available than the one performed here, the numerical examples 
indicate, that under still moderate rate constraints 1 — «; « 0.9 
holds, thus a = ^"l , w -^7= in (15). 

' \/2(l-«) 0.9\/2 ^ ' 
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